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INTRODUCTION 

In  response  to  a i^ques);?  from  Applied  Physics  Laboratory  at  Johns 
Hopkins  Unlversity'^or  a dynamic  analysis  of  an  underwater  vehicle  de- 
sign, NCSL  conducted  an  Investigation  of  the  Gibson's  Towed  Profiling 
System  (GTOPS)  vehicle. 

The  vehicle  consisted  of  two  like  bodies  Joined  by  a wing  and  tall 
appendages.  The  configuration  of  the  bodies  was  set,  but  the  size  and 
configuration  of  the  appendages  were  subject  to  change.  The  vehicle  was 
sized  and  analyzed  to  determine  longitudinal  and  lateral  stability. 

The  analysis  Included  the  effects  of  geometric  changes,  tow  cable 
considerations,  twin  fuselage  interference  drag,  the  effects  of  depth- 
keeping tethers,  horizontal  tall  actuator  sizing,  and  time  domain  dynamic 
vehicle  response.  Each  of  these  considerations  is  discussed. 

j'  'U.  ] 

VEHICLE  CONFIGURATION 


CONFIGURATION  CONTROL 

The  external  configuration  of  the  GTOPS  vehicle  had  been  controlled 
primarily  by  the  following  factors: 

1.  Payload  Characteristics.  Each  body  would  house  an  axisym- 
metrlc  Instrumentation  device.  The  two  bodies  are  Joined  by  a wing 
allowing  enough  separation  distance  to  avoid  interference  between  the 
instrumentation. 

2.  Variable  Depth/Towed  Operation  Requirement.  The  vehicle  must 
have  the  capability  to  change  or  hold  towing  depth  via  remote  control. 

3.  Dynamic  Stability.  Due  to  the  nature  of  the  measurements  to  be 
taken,  high  damping  ratios  and  low  modal  frequencies  were  desirable 
vehicle  characteristics. 

4.  Drag  Minimization.  In  order  to  permit  smooth  operation  of  the 
roller/cable  assembly,  the  vehicle  must  not  distort  the  cable  catenary 
in  the  region  of  the  towpoint. 
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Factors  1 and  4 led  to  the  twin  hull  GTOPS  configuration  using  a 
streamlined  form^^^.  Two  nearly  identical  body  forms  were  chosen  for 
analysis.  Both  forms  were  dependent  upon  the  size  and  shape  of  the 
instrumentation  packages  to  be  housed  in  the  bodies.  The  geometric 
characteristics  of  the  two  forms  were: 


Position  of 


Maximum 

Cross-Section 

X 

m 

Nose 

Radius 

R 

n 

Tail 

Radius 

R 

t 

Prismatic 

Coefficient 

C 

P 

Fineness 

Ratio 

L/D 

0.36 

0.7 

0.3 

0.65 

6.55 

- 655 

Body 

0.36 

1.0 

0.3 

0.65 

7.51 

- 751 

Body 

The  655  hull  shape  was  chosen  for  the  baseline  body  configuration 
in  the  Interest  of  reducing  size  and  weight. 


FORWARD  WING 

The  envelope  of  the  hardware  to  be  housed  Inside  of  the  forward 
wing  led  to  the  selection  of  an  0021  airfoil  section  shape  with  a chord 
length  of  1.5  feet.  The  distance  between  the  bodies  (root  chord  to  root 
chord  at  maximum  body  cross-sectional  area)  was  2.0  feet,  giving  a body 
centerline  separation  of  2.92  feet.  This  distance  was  chosen  with  re- 
gard to  the  lifting  force  required  of  the  forward  surface  for  depth 
change  response. 


ELEVATOR 

The  elevator  design  was  dependent  upon: 

1.  control  effectiveness  required, 

2.  longitudinal  stability,  and 

3.  hinge  moment /actuator  sizing. 

The  basic  design  criteria  was  to  keep  the  size  of  the  elevator  as 
small  as  possible  while  still  providing  adequate  effectiveness  and 
stability.  This  criteria  was  Intended  to  minimize  hinge  moments  and 
vehicle  unsteadiness  due  to  actuator  slop.  An  all  moveable  elevator  de- 
sign with  a chord  of  0.5  feet  and  a section  shape  of  an  0020  airfoil  was 
selected.  Shape  and  size  of  the  vertical  fins  were  almost  entirely 
dependent  upon  stability  considerations. 


^^^Davld  Taylor  Model  Basin  Report  719  Mathematical  Formulations  of 

Bodies  of  Revolution,  by  L.  Landwebcr  and  M.  Gertler,  September  1950. 
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FINAL  VEHICLE  DESIGN 

The  foregoing  brief  lilscusslon  of  the  GTOPS  vehicle  gives  a gen- 
eral Idea  of  the  design  criteria  used.  The  variational  analyses  pro- 
vide further  design  information.  The  resulting  vehicle  design  is  siiown 
in  Figures  I,  2,  and  3.  The  mass  distribution  data  are  given  in  Table  1. 


TABLE  1 

MASS  DISTRIBUTION  DATA 


Vehicle  displaced  weight  (lb) 

373.0 

3 

Vehicle  displaced  volume  (ft  ) 

5.82 

Body  wetted  area  (ft^) 

26.92 

Longitudinal  location  of  c.g.  (ft) 

(from  nose) 

2.t)6 

Longitudinal  location  of  c.b.  (ft) 

(from  nose) 

2.66 

Vertical  distance  c.b.  to  c.g  (ft) 
(positive  down) 

0.0417 

2 

X-axis  inertia  moment  (slug-ft“) 

9.0328 

2 

Y-axis  inertia  moment  (slug-ft*") 

7.4910 

2 

Z-axls  Inertia  moment  (slug-ft“) 

16.189 

GEOMETRIC  VARIATIONS  AND  STABILITY  ANALYSIS 


In  order  to  determine  the  effectiveness  of  the  base  case  design, 
several  computer  analyses  were  made  on  the  GTOPS  vehicle  with  variations 
from  the  base  case.  By  analyzing  the  stability  of  these  modified  vehicles 
a sensitivity  analysis  of  the  base  case  geometry  was  obtained.  T.ible  2 
describes  both  the  type  and  magnitude  of  base  case  variations  cv'nsiderevi. 


(Text  Continued  on  Page  71 
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roots  with  negative  real  parts).  Both  the  longitudinal  and  lateral 
roots  for  each  of  the  geometric  variations  are  given  in  Table  3.  Denomi- 
nator and  numerator  roots  for  each  of  the  longitudinal  and  lateral  trans- 
fer functions  are  presented  with  gain  values  in  Tables  4 and  5.  The 
stability  of  the  base  case  is  considered  along  with  a discussion  of  the 
effects  of  geometric  variations. 

BASE  CASE 


The  base  case  vehicle  proved  to  be  stable  in  the  longitudinal  domain 
and  neutrally  stable  in  the  lateral  domain.  The  longitudinal  roots  were 
all  real  indicating  that  no  oscillatory  modes  of  motion  were  present  (the 
base  considered  a cable  but  no  tethers).  The  existence  of  only  periodic 
roots  means  that  when  the  vehicle  is  disturbed  longitudinally  it  will  re- 
turn to  equilibrium  without  overshoot  or  oscillatory  motion.  The  lateral 
motion  was  also  non-oscillatory  and  stable  for  all  modes  except  roll  where 
it  was  only  neutrally  stable  due  to  a lateral  root  of  zero  (See  Table  5)  . 
The  lateral  roll  root  at  the  root  locus  origin  indicated  that  there  is 
no  mechanism  to  return  the  vehicle  to  a zero  roll  angle  vrfien  a lateral 
disturbance  is  encountered  due  to  the  zero  cb-cg  separation.  Sketches  of 
both  the  pitch  and  yaw  root  loci  are  shown  in  Figures  4 and  5 for  various 
cb-cg  separations. 

VARIATIONS 


Figure  4 shows  the  root  locus  plot  of  the  pitch  transfer  function  for 
changes  in  cb-cg  separation.  Three  cb-cg  locations  were  analyzed,  eg  V 
below  the  cb,  eg  at  the  cb,  and  eg  V'  above  the  cb.  The  pitch  transfer 
function  consists  of  four  poles  and  two  zeros^*^^^^.  Three  of  the  poles 
and  two  of  the  zeros  are  not  affected  by  cb-cg  position.  The  one  pole 
that  does  change  is  the  characteristic  mode  of  the  system  (l.e.  the  pole 
closest  to  the  origin);  therefore,  the  cb-cg  position  directly  affects  the 
response  time  of  the  vehicle.  For  the  eg  V below  the  cb  this  pole  is 
located  at  (s  + 0.057),  yielding  a settling  time  in  pitch  of  to.95  = 115 
seconds.  Locating  the  eg  V above  the  cb  moves  the  pole  into  the  right- 
hand  plane  to  (s  - 0.0045),  yielding  an  unstable  vehicle.  It  should  be 
noted  that  when  one  compares  the  time  histories  presented  later  in  the 
report  with  the  settling  time  of  115  seconds,  there  appears  to  be  a dis- 
crepancy; the  time  histories  show  very  rapid  response  to  pulse  and  doublet 


^’^Naval  Coastal  Systems  Laboratory  Report  287-76,  Development  of  the 
Equations  of  Motion  and  Transfer  Functions  for  Underwater  Vehicles, 
by  D.  E.  Humphreys,  July  1976. 

^^^Naval  Coastal  Systems  Laboratory  Report  242-76,  The  Analysis  of  a Longi- 
tudinal Control  System  for  Underwater  Vehicles,  by  D.  E.  Humphreys, 

R.  W.  Miller  and  L.  F.  Dewberry,  October  1975. 
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TABLE  4 

LONGITUDINAL  DENOMINATOR  AND  NUMERATOR  ROOTS 
WITH  TRANSFER  FUNCTION  GAINS 


Characteristic 

Equation 


Gain 

Roots 

-0.0260 

-0.2545 

-4.2294 

-11.486 

2.2754 

-0.0277 

-4.2558 

-10.724 

11.269 

-0.0995 

-0.2545 

-7.6360 

17.067 

-0.2545 

-3.9964 

11.269 

-0.2546 

-5.2091 

+7.6985 

TABLE  5 

LATERAL  DENOMINATOR  AND  NUMERATOR  ROOTS 
WITH  TRANSFER  FUNCTION  GAINS 


Characteristic** 

Equation 


6.0237 

7.0972 

0.0 

-6.0237 


*When  the  base  case  data  were  modified  to  produce  non-zero  roll  numerator 
coefficients  It  was  found  that  the  roll  transfer  function  contained  a 
free  zero  root  In  the  denominator. 

**A  5th  root  exists  In  the  lateral  CE  with  a value  of  zero  which  exactly 
cancels  with  a free  zero  In  each  lateral  transfer  function. 
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FIGURE  U.  ROOT  LOCUS  DIAGRAM  OF  VEHICLE  PITCH  FOR  CB-CG  VARIATION 
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inputs.  The  rapid  response  shovm  in  the  time  histories  is  due  to  the 
selection  of  an  all  movable  horizontal  tall  rather  than  a fixed  plane 
with  only  a small  portion  movable.  This  results  in  a high  gain  system. 
For  example,  a 5-degree  step  in  the  horizontal  tail  produces  a steady 
state  pitch  angle  of  119  degrees  with  a settling  time  of  115  seconds. 

For  this  same  input  after  only  5.8  seconds,  the  pitch  angle  would  be 
28.6  degrees.  Thus,  large  pitch  motions  can  be  achieved  rapidly  with 
small  tall  deflections. 

The  effect  of  cb-cg  variation  on  the  lateral  dynamics  is  shown  in 
the  yaw  root  locus  of  Figure  5.  The  yaw  transfer  function  consists  of 
three  zeros  and  five  poles.  Only  two  of  the  poles  and  two  of  the  zeros 
are  affected  by  cb-cg  change.  For  the  eg  V below  the  cb  the  poles  are 
oscillatory  with  a damping  ratio  of  0.12  and  an  undamped  natural  fre- 
quency of  0.36  radians /second.  With  the  cb  at  the  eg,  the  poles  are 
located  at  (s  + 0.0)  and  (s  + 0.089).  For  the  eg  V above  the  cb,  the 
poles  are  at  (s  + 0.41)  and  (s  - 0.32)  indicating  an  unstable  system. 

The  two  zeros  of  the  yaw  transfer  function  that  are  changed  by  cb-cg 
variation  are  seen  to  closely  track  the  movement  of  the  poles,  effec- 
tively canceling  the  contribution  of  these  poles  on  yaw  response.  Be- 
cause of  this  cancellation,  the  net  effect  of  cb-cg  variation  on  yaw 
response  is  zero.  This  same  result  was  observed  in  the  side  velocity 
transfer  function.  However  for  the  roll  transfer  function,  the  zeros 
do  not  cancel  the  poles,  and  there  is  a pronounced  effect  of  cb-cg  vari- 
ation on  roll  response  as  shown  by  the  lateral  time  histories  presented 
later  in  the  report. 


Varying  the  distance  between  the  nose  and  aerodynamic  center  of  the 
wing  produced  little  movement  of  the  stability  roots  with  the  tow  point 
considered  to  remain  at  the  wing  aerodynamic  center.  The  aspect  ratio 
variation  of  the  horizontal  tail  also  had  little  effect  on  root  location 
indicating  little  change  in  overall  vehicle  dynamics.  Variation  in  the 
cb-cg  did  however,  produce  a significant  effect.  With  the  eg  above  the 
cb  the  vehicle  became  unstable  in  both  the  longitudinal  and  lateral  modes. 
The  low  center  of  gravity  configuration  introduced  an  oscillatory  roll 
mode  in  the  lateral  domain  but  removed  the  neutral  stability  of  the  base 
case  by  driving  the  lateral  root  at  the  origin  into  the  left-half  plane. 
Since  this  variation  proved  stabilizing  for  both  longitudinal  and  lateral 
dynamics  it  definitely  represents  an  improvement  to  the  base  case 
configuration. 

Little  variation  in  stability  was  noted  for  the  body-body  separa- 
tion, the  longitudinal  eg  position,  and  the  weight  and  buoyancy  variation, 
as  noted  by  Table  3 which  shows  the  lateral  and  longitudinal  roots.  The 


(Text  Continued  on  Page  14) 
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FIGURE  5.  ROOT  LOCUS  DIAGRAM  OF  VEHICLE  YAW  FOR  CB-CG  VARIATION 
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vertical  tall  size  variation  did,  however,  produce  an  unstable  lateral 
mode  for  the  vertical  tall  size  which  was  50  percent  of  the  original 
area.  The  larger  tall  area  produced  a more  stable  vehicle  as  expected. 
The  size  of  the  base  case  vertical  tall  Is  sufficient  from  a stability 
point  of  view  but  should  not  be  reduced  by  a great  extent  due  to  the 
possible  unstable  mode. 

The  speed  variation  Indicated  that  the  vehicle  was  more  stable  at 
higher  speeds. 

The  variation  obtained  with  the  751  body  produced  little  or  no 
change  In  vehicle  dynamics. 


TWIN  FUSELAGE  INTERFERENCE  DRAG 


The  possibility  of  an  optimum  separation  distance  between  the  twin 
fuselages  for  minimum  drag  has  been  considered.  As  the  separation  dis- 
tance decreases,  the  Interference  drag  between  the  two  fuselages  in- 
creases, but  the  wing  drag  decreases  (since  Its  planform  area  decreases). 
Hence,  there  could  be  an  optimum  separation  distance  for  minimum  drag 
which  would  depend  on  the  manner  in  which  the  fuselages'  interference 
drag  varies  with  separation  distance.  Consider  the  sketch  shown. 


The  Interference  drag  is  a function  of  Y/d  and  is  obtained  from 
Hoerner's  "Fluid  Dynamic  Drag,"  page  13-16^^^.  Defining 

C ■ total  drag  including  Interference 
°T0T 

Cj^  ■ drag  at  separation  Y * infinity, 

INF 


(4) 


Hoerner,  S.  F. , Fluid  Dynamic  Drag,  Published  by  author,  1965. 
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the  variation  of  C_  /C_  versus  Y/d  Is  shown  in  Figure  6. 

TOT  INF 

Obviously,  the  optimum  separation  distance  will  be  at  some  value  of  Y/d 
less  than  0.5.  On  the  other  hand,  the  present  GTOPS  configuration  has 
a Y/d  ■ 1.09.  If  the  separation  distance  were  reduced  to  Y/d  less  than 
0.5,  the  loss  of  wing  span  (and  planform  area)  would  result  in  Insuffi- 
cient vehicle  lift,  unless  outboard  wings  were  added  which  would  in 
turn  Increase  the  drag.  Since  the  theoretical  optimum  design  is  not 
very  practical,  a separation  distance  was  chosen  which  would  yield  a 
reasonable  amount  of  wing  lift  rather  than  an  optimum  drag  conf igurat ion. 


FIGURE  6.  INTERFERENCE  DRAG 


i 

I 

; EFFECT  OF  DEPTH- KEEPING  TETHERS 

! 

ASSUMPTIONS 

i 

i Reference  5 has  been  reviewed  and  appears  to  be  correct  within  the 

I limitations  of  the  main  assumptions:  (1)  neutrally  buoyant  tethers,  (21 

; tangential  cable  drag  negligible,  and  (3)  main  tow  cable  straight  and 

vertical.  Therefore,  the  computed  effective  vertical  spring  constants 
have  been  used  In  this  analysis. 

In  this  analysis,  the  depth-keeping  tethers  are  modeled  as  a vertl- 
I cal  spring  of  stiffness  k attached  to  the  vehicle.  The  effect  of  this 

vertical  spring  appears  in  the  vehicle's  equations  of  motion  as  the  non- 
dimensional  hydrodynamic  coefficients  which  follow: 


^^^Applled  Physics  Laboratory,  The  Johns  Hopkins  University  Report  DM- 
3591,  GTOPS  .'-lechanical  Developwf>nt  Pi'ogress  Report  and  Program 
Documentation,  by  W.  A.  Vcnzla,  20  June  1977. 
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where 


Z'  - ■ 

2'  - Z’ 

* H PU^t 

6 z 

«9  - 

(X^p/O 

(X^p/O^ 


U - vehicle  speed 

i - reference  length  (-  6.0  ft  for  APL  body) 


longitudinal  distance  between  eg  and  attachment 
point  of  vertical  spring  (positive  ahead  of  eg) 


p ■ mass  density  of  sea  water. 


Table  6 shows  spring  constants  of  Reference  3 converted  to  values 
of  as  a function  of  tether  length  and  vehicle  speed.  Note  that  Z*  is 
essentially  independent  of  speed,  especially  for  the  longer  tether.  ^ 


TABLE  6 

SPRING  CONSTANTS  CONVERTED  TO  Z'  VALUES 

z 


Te  ther 
Length 
(ft) 

U 

(knots) 

k 

(Ibf  1 ft) 

Z'  X 10^ 
z 

64.2 

4 

1.8 

-1.884 

64.2 

8 

6.4 

-1.675 

64.2 

10 

10.4 

-1.742 

96.3 

4 

0.95 

-0.994 

96.3 

8 

3.75 

-0.981 

96.3 

10 

5.6 

-0.938 
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ANALYSIS 

Spring  at  CG,  No  CG-GB  Separation 

I Three  cases  were  analyzed.  The  first  considered  the  effect  of 

varying  the  spring  constant  k for  a vehicle  with  the  spring  attached  at 
I ^ the  eg  with  zero  cb-cg  separation.  The  vehicle  speed  was  4 knots  and 

' the  spring  constant  values  were  0.0,  0.4,  0,8,  and  1.6  Ibf/ft,  An 

I investigation  of  the  roots  indicated  that  a non-zero  spring  constant 

I introduces  an  oscillatory  longitudinal  mode  with  a decreasing  period  as 

k increases. 

‘ i 

j Spring  at  CG,  CG  Below  CB  j 

I The  second  case  covered  the  effect  of  varying  tlie  spring  constant  k 

t for  a vehicle  with  the  spring  attached  at  the  eg  and  the  eg  one-half  inch 

1 below  the  cb.  Again,  the  speed  was  4 knots  and  k had  the  values  0.0,  0.4, 

‘ 0.8,  and  1.6  Ibf/ft.  The  results  indicated  that  again  as  k increases 

' above  zero,  oscillatory  motion  was  encountered  in  the  same  manner  as 

< above;  however,  the  overall  damping  was  increased  due  to  the  effect  of 

cb-cg  Separation. 

' Variation  in  Tether  Attachment  Point  Location 

I The  third  case  analyzed  the  effect  of  tether  attachment  location. 

The  parameters  for  this  study  were  a speed  of  4 knots,  zero  cb-cg 
' separation,  a k of  1.6  Ibf/ft,  and  tether  attachment  point  locations 

(X^p)  of  1.0,  0.5,  0.25,  0.0,  -0.25,  -0.5,  and  -1.0  feet.  The  resulting 
roots  (Table  7)  indicate  that  the  vehicle  becomes  unstable  if  the  tether 
attachment  point  is  behind  the  eg  (negative  values  of  X.pp).  Figure  7 
j shows  the  effect  of  the  longitudinal  location  of  the  tethers  (X-pp)  on 

stability.  The  pitch  transfer  function  has  five  poles  and  three  zeros; 
note  that  an  additional  pole  and  zero  are  present  due  to  the  tethers. 

Only  two  of  the  poles  and  zeros  are  affected  by  X-pp  variation.  For 
0 ^ X'pp  ^ 1.0,  the  system  exhibits  lightly  damped  characteristics  which 
means  the  system  will  have  faster  response  and  more  overshoot  tlian  the 
system  without  tethers  due  to  the  spring  effect  of  the  tethers.  For  ex- 
ample at  Xpp  » 0.5,  the  poles  are  at  -0.019  ± j 0.12  and  the  zero  is  at 
-0.0166.  The  two  poles  are  the  characteristic  mode  of  the  system  since 
they  are  closer  to  the  origin  than  all  other  system  poles.  The  damping 
ratio  for  these  poles  is  0.16.  Because  the  zero  is  closer  to  the  origin 
than  the  poles,  the  system  will  exhibit  lead  characteristics.  For 
-1.0  £ Xpp  < 0,  l.e.  the  attachment  point  behind  the  eg,  the  system  is 
unstable,  as  witnessed  by  the  single  pole  in  the  right-half  plane. 


(Text  Continued  on  Page  19) 
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TABLE  7 

LONGITUDINAL  ROOTS  VOR  VARIOUS  TETHER 
ATTACHMENT  POINT  LOCATIONS 


Longitudinal  Denominator  Roots 


l.O 

-0.2546 

-4.2152 

-11.486 

(-0.020  i j 0.1587) 

0.5 

-0.2546 

-4.2169 

-11.486 

(-0.019  ♦ j 0.1213) 

0.25 

-0.2546 

-4.2177 

-11.486 

(-0.019  ± j 0.0973) 

0.0 

-0.2546 

-4.2185 

-11.485 

(-0.0188  t j 0.0650) 

0.25 

-0.2546 

-4.2190 

-11.484 

-0.0506 

+0.0131 

0.5 

-0.2546 

-4.2200 

-11.484 

-0.0976 

+0.0607 

1.0 

-0.2546 

-4.2220 

-11.485 

-0.1468 

+0.1121 

HORIZONTAL  TAIL  ACTUATOR 


The  torque  required  by  the  horizontal  tail  actuator  ir,  dependent 
upon  the  airfoil  characteristics  and  the  position  of  the  scock.  Posi- 
tioning the  stock  as  close  as  possible  to  the  center  of  pressure  will 
minimize  the  hinge  moment  due  to  the  normal  force  acting  at  the  center 
of  pressure.  However,  should  the  center  of  pressure  move  behind  the 
stock  a reversal  in  the  hinge  moment  will  occur.  Depending  on  the  hold- 
ing torque  characteristics  of  the  actuator,  this  phenomenon  can  cause 
the  control  surface  to  rotate  to  an  uninitiated  extreme  deflection 
angle.  The  stock  location  is  generally  placed  somewhat  aft  of  the  pre- 
dicted locus  of  the  center  of  pressure  to  avoid  this  problem.  For  this 
design  the  location  of  the  center  of  pressure  is  expected  to  vary  between 
20  and  23  percent  of  the  chord  length  for  nonstall  angles  of  attack 
(|a|  ^ 25  degrees).  A stock  position  at  the  25  percent  chord  yields 
hinge  torques  required  to  balance  the  moment  due  to  the  normal  force 
of  between  1600  and  1800  inch-ounces  for  a = 10  degrees  and  a speed  of 
12  knots.  At  4 knots  this  reduces  to  between  180  and  200  inch-ounces. 

It  should  be  recognized  that  an  angle  of  attack  of  10  degrees  is  quite 
moderate.  Required  hinge  torques  will  go  up  somewhat  less  than  directly 
proportional  to  the  angle  of  attack. 

An  important  consideration  uere  is  that  the  holding  torque  must  be 
sufficient  to  keep  the  control  surface  at  a desired  angle  or  vehicle 
characteristics  will  be  altered. 
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KONCITUDINAL  TIME  DOMAIN  ANALYSIS 


TEST  CONFIGURATIONS 


Tilt’  effects  of  cb-cg  separation  and  depth-keeping  tethers  on  the 
vehicle's  longitudinal  dynamics  were  determined  from  time  histories. 

These  time  histories  were  generated  by  exciting  the  longitudinal  trans- 
fer functions  of  the  various  configurations  with  various  elevator  forcing 
functions.  Samples  of  the  elevator  or  stern  plane  forcing  functions  are 
presented  In  Figure  8.  The  longitudinal  variables  examined  In  this  se- 
ries of  tests  were  axial  velocity  u,  normal  velocity  w,  pitch  angle  i', 
and  depth  change  z.  The  matrix  of  tests  as  characterized  by  tether 
S' Iffness,  cb-cg  separation,  and  elevator  forcing  function  Is  shown  In 
Tai.l-'  8. 

TABLE  8 

LONGITUDINAL  DYNAMICS  TESTS 


Elevator  Function 

5 deg, 

1 sec  Pulse 

5 deg, 

5.8  sec  Pulse 

5 deg, 

0.5  sec  Doublet 

5 deg. 

1 sec  Pulse 

5 deg. 

5.8  sec  Pulse 

5 deg. 

0.5  sec  Doublet 

5 deg. 

1 sec  Pulse 

5 deg, 

5.8  sec  Pulse 

3 deg. 

0.5  sec  Doublet 

5 deg. 

1 sec  Pulse 

5 deg. 

5.8  sec  Pulse 

5 deg. 

0.5  sec  Doublet 

15  deg. 

1 sec  Pulse 

15  deg. 

5.8  sec  Pulse 

15  deg. 

0.5  sec  Doublet 
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RKSULTS 


TVo  polntN  MhiMit  the  roMultH  urovi  to  I'o  ompli.-is  1 . I’ltst,  t l\i- 

roMultii  |tlv9i\  h«*rt'  jiro  for  a volilolo  npoiHl  ol  •*  Roots.  I'ho  t t on, Is  sluoiM 
bo  the  aumo  «t  otbor  npooilH,  bol  t bo  m.innt  t islos ; o.y\.,  .loi>tb  ob.tuno.  will 
i\ot . Sov'Oiul,  Momi>  oiiMoM  liwolvlt\>i  pulsos  v)t  soooiuls  Jvn.it  li'ii  .ii  v'  In- 

valid bocAUH*  thov  tnvolvo  vv'bloU'  iiu't  iv»ns  ti'spoi  1 .1 1 I v pltilO  wblob  osvv'wl 
Ibo  amall  portnrb.K  ton  assumptions  >>11  wlitoli  tbo  llno.ti  .in.tivsis  Is  b.isoil. 

W 1 1 h o u t jr  a_t  h o r a 

Wlthv>ut  tothovs,  tbo  V4'btv'l«>  iiu't  li'n  Is  ossv'nt  1 .1 1 I v apoi  1 v'vl  K- , .nnl 
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Ib'wovor,  oh-ojt  Moparat  IvMi  (i'k  b«>li'w  i lO  Im'i'o.'isos  tbo  plii'li  sl.ibllliv 
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oxamplo,  tho  depth  v'han)ti'  oansod  bv  .1  '»  d«'n,  1 siS'  I'vvlso  Is  lovivn'vsl  I i om 
■lA  foot  to  .*8  loot  wbon  tin'  o>t  Is  iiiovod  1/.’  Im  b bi'low  I lu'  i b. 

With  Tothora 


With  tothoi's,  t b«'  vv'blolo  iiuvl  K'n  Is  ivsv- 1 1 l.it  I'lv  Iwltb  tbo  oxv'0)'t  Iv'n 
ol  n whioh  romalns  .ipoi  It'il K'  b«'i'anso  tbo  vi'i  tU'.'il  spi  ln>t  ib'os  not  v'onpio 
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Kv>r  a tv’thor  atlltnosa  ol  I . l<  Ibl/tl  I ' > , lb*'  il.impod  n,iliii';il 

porliul  la  Itll.S  Hoovvnd  aiul  Ibo  I liiu'  to  damp  Is  88  sov  v’nd  t'  I Imo  oon- 
stanta),  wbllo  tv'i  a tolbor  sllltnv'sii  v'l  !'<  Ibl/ll  0'’l>tin<'s  1.'  .nivi  1 O, 
tho  damped  natural  porliul  Is  U./  soi-oiul  anil  tin'  I Imo  1 1>  iljimp  Is  'ts 
aooond.  It  abonid  ho  m'tovl  that  Ibo  slittnoss  .at  s Riu'Is  will  |irv>bablv 
bo  on  tho  I'rdor  I't  l.n  Ibf/lt  ami  m'l  I ■'  Ibl'll. 
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IK.ST  I'.ONFl DURATIONS 
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torlaad  by  ob-cji  aoparat  Ion,  rndilor  lvp«',  .nul  nnbloi  li'ii  lnn  Ivnn'ili'ii  Is 
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and  U'wer  vortical  atabl  I Icors.  This  spovlal  oas«'  was  inn  simpiv  tv'  pro- 
vide a moana  of  exciting  tbo  roro  I'b-v  ^ sop<irat  h'li  I'oiil  t^iii  .it  Iv'ii  In  .>'11. 
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FIGURE  9.  LONGITUDINAL  TLME  HISTORIES  FOR  TEST 
NTMEERS  1,  2,  AND  3 OF  TABLE  8 


FIGURE  10.  LONGITUDIHAL  TIME  HISTORIES  FOR  TEST 
NUMBERS  4,  5,  AND  6 OF  TABLE  8 
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TABLE  9 

LATERAL  DYNAMIC  TESTS 


Test  No. 

cb-cb  Separation 
(in.) 

Rudder  Type 

Rudder  Function 

1 

0.0 

No  Lower 

5 deg,  1 sec  Pulse 

2 

0.0 

No  Lower 

5 deg,  0.5  sec  Doublet 

3 

0.0 

Symmetric 

5 deg,  1 sec  Pulse 

4 

0.0 

Symmetric 

5 deg,  0.5  sec  Doublet 

5 

0.5 

Symmetric 

5 deg,  1 sec  Pulse 

6 

0.5 

Symmetric 

5 deg,  0.5  sec  Doublet 

RESULTS 


As  noted  in  the  longitudinal  analysis,  the  results  given  here  are 
for  4 knots.  The  trends  should  be  the  same  at  other  speeds,  but  the 
magnitudes;  e.g.,  side  track,  will  not. 


This  analysis  assumes  that  neither  the  tow  cable  nor  the  depth- 
keeping  tethers  produce  roll  moments  or  side  forces.  In  the  actual 
vehicle,  they  will  be  coupled  with  both  roll  and  side  track,  but  deter- 
mination of  these  coupling  terms  would  require  a more  sophisticated 
analysis.  The  present  analysis,  however,  has  demonstrated  the  need  to 
have  both  roll  and  side  track  restoring  moments  and  forces  present  in 
adequate  amounts.  The  time  history  plots  described  in  Table  9 are  pre- 
sented in  Figures  14  through  16. 


No  cb-cg  Separation 


Without  cb-cg  separation,  the  vehicle  motion  is  essentially 
aperiodic  and  any  disturbance  will  produce  finite,  steady  state  values 
of  side  track  and  roll  (Figure  14).  For  example,  a 5 degree,  1 second 
pulse  of  only  the  upper  rudder  causes  a final  roll  angle  of  16  degrees 
and  a final  side  track  of  2.4  feet. 
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FIGURE  14.  LATERAL  TIME  HISTORIES  FOR  TEST  NUMBERS  1 AND  2 OF  TABLE  9 


STORIES  FOR  TEST  NUMBERS  3 SNF  4 OF  TABLE  9 
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With  cb-CK  Separation 

With  cb-cg  separation,  the  vehicle's  roll  motion  becomes  oscilla- 
tory (v,  lii,  and  y remain  aperiodic)  and  the  final,  steady  state  angle 
will  be  zero  (Figures  15  and  lb).  However,  there  will  still  be  a finite, 
steady  state  side  track.  For  1/2-lnch  cb-cg  separation,  the  damped 
natural  roll  period  Is  17.4  seconds  and  the  time  to  damp  is  68  seconds 
(3  time  constants). 


CONCLUSIONS 


1.  The  mass  distribution  of  the  vehicle  should  be  designed  so  that 
the  center  of  gravity  is  positioned  vertically  below  the  center  of  buoy- 
ancy as  far  as  possible.  This  will  provide  roll  and  pitch  stability  as 
well  as  minimum  depth  loss  due  to  disturbances. 

2.  A horizontal  tail  actuator  capable  of  delivering  50  inch-ounces 
of  torque  is  insufficient  for  moderate  elevator  deflection  angles.  At  a 
speed  of  4 knots  approximately  100  inch-ounces  will  be  required  for  a 
deflection  angle  of  5 degrees  (this  is  for  the  tethered  vehicle  with 
1/2  in.  cb-cg  separation).  An  actuator  slop  of  1/4  of  a degree  for  the 
elevator  will  result  in  an  error  of  1 degree  in  pitch  and  8 feet  in  depth 
at  4 knots. 

3.  The  stock  position  for  the  elevator  should  be  located  no  further 
forward  than  the  25  percent  chord  position  of  the  elevator. 

4.  The  tethers  tend  to  drive  the  depth  changes  to  zero  for  longi- 
tudinal disturbances.  An  elevator  control  system  would  be  required  for 
no  depth  loss  without  tethers. 

5.  The  tethers  should  be  attached  ahead  of  the  vehicle  eg  because 
their  attachment  aft  of  the  eg  would  cause  the  vehicle  to  be  unstable. 

6.  Our  simplified  model  predicts  only  longitudinal  effects  for  the 
tethers;  the  lateral  effects  could  be  estimated  with  a more  in-depth 
analysis.  They  do  appear,  however,  to  provide  a stabilizing  effect  in 
roll  and  yaw. 

7.  The  tether-vehicle  analysis  assumes  that  there  is  no  signifi- 
cant interaction  between  the  tethers  and  the  depressor  below  the  body 
(the  depressor  is  assumed  to  institute  no  dynamic  disturbances  to  the 
system) . 
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8.  Based  on  the  tether  analysis  provided  to  NCSL^‘’\  a 1.6  pound/ 
foot  spring  constant  is  a reasonable  estimate  at  4 knots.  The  frequency 
of  oscillation  predicted  by  this  analysis  is,  however,  too  high  by  a 
factor  of  the  square  root  of  32.17.  The  analysis  presented  in  this  report 
for  a spring  constant  of  1.6  calculated  a longitudinal  oscillatory  period 
of  101.5  seconds.  This  period  is  much  higher  than  that  predicted  by  the 
analysis  treating  the  body  as  a spring  mass  system  assuming  no  pitch 
coupling.  This  difference  implies  that  the  coepling  cannot  be  neglected. 
In  general  it  is  questionable  whether  or  not  a linear  spring  constant  of 
15  pound/foot  can  be  achieved  for  the  anticipated  tethered  vehicle  design. 

9.  Although  a vertical  cb-cg  separation  does  provide  for  roll  sta- 
bilization, some  consideration  should  be  given  to  adding  partial  span 
ailerons  to  the  wing  of  the  body  to  provide  adequate  roll  control. 

10.  Coupling  of  longitudinal  and  lateral  motions  are  instituted  by 
both  hydrodynamic  and  tether  forces.  A nonlinear  analysis  would  require 
the  knowledge  of  nonlinear  hydrodynamic  and  coupling  coefficients.  The 
linear  analysis  provided  here  assumes  no  coupling.  A nonlinear  analysis 
would  require  a much  more  detailed  Investigation.  Past  experience  with 
other  towed  vehicles  has  Indicated  that  as  long  as  vehicle  motions  about 
an  equilibrium  state  are  "small"  the  linear  equations  have  been  found  to 
adequately  describe  its  dynamics. 


